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ABSTRACT 

Poynting-Robertson drag has been considered an ineffective mechanism for delivering 
dust to regions interior to the cool Kuiper belt analogues seen around other Sun-like 
stars. This conclusion is however based on the very large contrast in dust optical depth 
between the parent belt and the interior regions that results from the dominance of col¬ 
lisions over drag in systems with detectable cool belts. Here, we show that the levels 
of habitable zone dust arising from detectable Kuiper belt analogues can be tens to a 
few hundreds of times greater than the optical depth in the Solar Zodiacal cloud. Dust 
enhancements of more than a few tens of ‘zodi’ are expected to hinder future Earth¬ 
imaging missions, but relatively few undetectable Kuiper belts result in such levels, 
particularly around stars older than a few Gyr. Thus, current mid to far-IR photomet¬ 
ric surveys have already identified most of the 20-25% of nearby stars where P-R drag 
from outer belts could seriously impact Earth-imaging. The LBTI should easily detect 
such warm dust around many nearby stars with outer belts, and will provide insight 
into currently unclear details of the competition between P-R drag and collisions. Given 
sufficient confidence in future models, the inevitability of P-R drag means that the non¬ 
detection of warm dust where detectable levels were expected could be used to infer 
additional dust removal process, the most likely being the presence of intervening plan¬ 
ets. 

Key words: circumstellar matter — zodiacal dust — planets and satellites: detection 
— radiation: dynamics 


1 INTRODUCTION 

As the regions deemed most likely to harbour alien life, the hab¬ 
itable zones of other stars are becoming of increasing interest. 
Major efforts to find p otentially habitable p l anets are underwa y 
(e.g. Kepler , PLATO. |Borucki et alJl2003l : iRauer et al.1120141) . 
with the goal of pushing the detection limits towards true Earth- 
analogues. Coupled to this interest has been parallel work on 
exo-Zodiacal dust, both as a potentially useful dynamical tracer 
that reveals unseen planets (e.g. Sta rk & KuchneJl2008h . and as 
a source of noise and confusion that may hinder the direct de- 
tection and characterisation of such planets (Defrere et alj20icl : 
iRoberge et al.ll2012l : IStark et aill20 1 4l : lBrownll2015h . 

Habitable zone (HZ) dust is indeed seen around other stars 
(e.g. ISong et al .1 lioosl : iFuiiwara et alj|2010h . though these de¬ 
tections are limited by the photo metric methods used to r ela- 
tively high dust levels that are rare dKennedv & Wvattl2013l ). In 
general the origin of this exo-Zodiacal (‘exo-zodi’) dust is un- 
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known. In some cases, particularly for main-sequence stars, the 
levels are sufficiently extreme that the dust cannot simply orig¬ 
inate in massive analogues of our Asteroid belt because colli- 
sional decay would have ground the d ust le vels well below those 
observed over the stellar lifetime dWvatt et~afl l2007ah . Possi¬ 
ble solutions are th at the dust was created in a recent c ollision 
(i.e. is transient, e.g. ISong et al .120051:1 Mena et al.l2012l) . or that 
the dust is being continuously delivered in the form of comets 
from more distant regions (e.g. (Wyatt et al. 2007a; B onsor et alJ 
120121 ) . 

Curiously, the mechanism responsible for delivering much 
of the Zodiacal dust from exterior regions to the Earth’s vicin¬ 
ity, Poynting-Robertson (P-R) drag, can be ruled out in almost 
all extra-Solar cases. The reason is that particles spiralling in by 
P-R drag are also subject to collisional destruction and subse¬ 
quent radiat ion-pressure blowout (for stars more luminous than 
the Sun, see lReidemeister et al.il201 ll) . The fate of dust created 
in a parent belt of asteroids or comets (i.e. a ‘debris disk’) there¬ 
fore depends on the relative importance of P-R drag and colli¬ 
sions. This issue was explored by Iwvatt ( 2005 ). who used a 
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simple model_of P-R drag for dust of a single particle size (see 
also lWvatt et alJfl999l ). The key points were: i) the denser the 
parent belt, the greater the contrast in optical depth between 
the parent belt and the inner regions, ii) the denser the parent 
belt, the greater the absolute dust level in the inner regions up to 
some maximum level, and iii) the debris disks that are detected 
around other stars are generally collision dominated, rather than 
P-R drag dominated, and hence P-R drag could be considered 
insignificant. The motivation was to show that the invocation of 
extra forces (e.g. unseen planets) was not necessary to explain 
why the regions interior to some cool disks were not filled in. 
That is, the model was largely considered within the context of 
the relative levels of dust in the parent belts and in the interior 
regions. 

_A n aspe ct of the model that was not explored in detail by 

Iwvatlj 12003 ) was the absolute level of interior dust, in part be¬ 
cause this was not the primary motivation, but also because ob¬ 
servations at the time were not sensitive enough to detect the 
dust levels predicted. However, new instruments are pushing the 
sensitivity limits for warm dust to levels low enough that P-R 
drag as an origin of habitable zone dust is becoming relevant. 
While the photometric methods used to discover bright dust are 
limited to mid-IR disk to star flux ratios greater than about 10% 


(at 3c r), the Keck Interferometer Nuller (KIN, Serabyn et al 


2012) was able to reach levels closer to 1% |\4 i 11an-Gab et et al 
201 ll : iMennesson et alj|2014lt, and the Large Binocular Te le- 


scope Interferometer (LBTI, Hinj|200^ : lDefrere et alJl2015h is 

expected to go be low 0.1%. _ _ 

Specifically, iMennesson et alj l i2014l) report new results 
from the KIN, showing that the level of warm dust seen around 
nearby stars tends to be higher in systems with detections of 
cool Kuiper belt analogues. The surprising aspect is that in most 
cases the levels tend to be similar, a few hundred times the Solar 
System level. The typical explanation for a direct link between 
warm and cool belts is comet delivery by p lanet scattering (e.g. 
iBonsor & Wvattll20l2l : lBonsor et alJl2012l) . but the outcome of 
this scenario depends sensitively on various parameters such as 
the planetary system architecture, so is an unlikely explanation. 
In contrast, the competition between P-R drag and collisions 
naturally leads to warm dust levels that are relatively insen sitive 
to the properties of the parent belt, and lMennesson et alj ( 1201 4) 
show that their results are consistent with levels expected from 
the lWvatll d2005l) model of P-R drag delivery 

Here, we consider t he wider implicatio n s of th ese findings, 
including the results of Ivan Lieshout et alj (120141) who devel¬ 
oped a more sophisticated numerical model of the competition 
between P-R drag and collisions. Their res ults a re easily incor¬ 
porated into the analytic model of lWvattl J2005I ). which we use 
to make general predictions of the warm dust levels in systems 
with Kuiper belt analogues. 


2 MODEL 

2.1 Optical depth interior to parent belt 

The steady-state P-R drag model of lWvatll l2005ll calculates the 
dust optical depth as a function of radius interior to a source 
region, a parent belt of planetesimals where dust is created in 
collisions between larger bodies. It is an analytic solution of 
the continuity equation, where the number of particles enter¬ 
ing a radial region equals that leaving plus that lost to colli¬ 


sions. The main assumptions were that particles have a single 
size, and that all collisions were destructive and result in per¬ 
manent loss of the fragments due to stellar radiation pressure. 
The model does not therefore apply to stars with insufficient lu¬ 
minosity or stellar wind pressure to blow small dust out of the 


system (see Augereau & Beust 2006; 

Reidemeister et al. 2011; 

Lestrade et al. 20l3 Schiippler et al. 

2014], for discussions of 


this case). 


Particles at some outer radius r 0 , representing the parent 
belt, have a face-on geometric optical depth to, and the only re¬ 
maining free parameters are the stellar mass and the strength of 
the P-R drag on the particles. The latter is parameterised by the 
ratio of the force felt by a particle due to radiation pressure to the 
gravitational force, /3 = F ra d/ F* oc L*/M* (i.e. we assume 
that the analogous stellar wind drag effect is relatively small). 
Small particles created on initially circular orbits are blown out 
of the system on a dynamical timescale when /3 > 0.5, so the 
model assumed /3 = 0.5. 

A numerical versio n of this model was developed by 
ivan Lieshout et all 120141) . who included a dust size distribution 
and accounted for the fact that grains of different sizes are af¬ 
fected to differing degrees by both radiation pressure and P-R 
drag, and also have different collisional lifetimes due to differ¬ 
ences in orbital eccentricity, strength and their relative numbers. 
Compared to the analytic model, the only significant difference 
was that the dust levels in the regions interior to the parent belt 
were nearly an order of magnitude lower. The reason is that 
the smallest bound grains are created on highly eccentric orbits 
meaning that collisions are more frequent and at higher rela¬ 
tive velocities, so a given particle can be destroyed by a smaller 
(more common) impactor. They found that interior to the par¬ 
ent belt, the size distribution is dominated by grains near the 
blowout size, which is consistent with the single grain size as¬ 
sumed in the analytic model. They only explored cases in which 
the parent belt was relatively massive, so it seems likely that 
this conclusion arises due to the dominance of collisions imme¬ 
diately interior to the parent belt, where larger grains are heav¬ 
ily depleted by the more numerous near-blowout size grains. In 
the case where collisions are much less important larger grains 
would be expected to dominate because of their slower P-R drag 
timescale. Our focus here is to relate observations to detectable 
parent belts, so the approximation of a single grain size at the 
blowout limit is reasonable, but more numerical work is needed 
to see how the results change with the parent belt optical depth. 

For relatively dense parent belts a simple way to include 
the effect of elliptical orbits in the analytic model is to change 
the collision timescale by a multiplicative factor, k. The solution 
to the continuity equation for the optical depth, r(r), is 

r(r) = r 0 / [l + 4 770 (l - v/r/ro^j , ( 1 ) 

rjo = 5000T O ^(r o /ae)(M 0 /M*)/(/3fe). (2) 


These equations only differ from those inlWvatt 1 2005|) by the 
factor k. Given the results of lvan Lieshout et alj 20141) , we can 
assume that /3 = 0.5 because these grains domina te th e si ze dis - 
tribution. We set k = 1 to represent the original Iwvatll (1200 50 
model, which is consistent with the KIN results and is referred 
to as the “low co llision rate”, or set k = 1 /7 to match the nu¬ 
merical model of lvan Lieshout et~afl ( l2014h . which we refer to 
as the “high collision rate”. In the absence of the KIN results 
we would prefer the high collision rate because it was derived 
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by a model thought to be more realistic. However, an apparently 
lower collision rate could arise because the numerical model 
does not include all important physical effects, for example if 
not all collis ions result in destruction and subsequent blowout 
(e.g. lKriit & Kamall2014l) . One of our goals is therefore to make 
predictions for each case so that observations can empirically 
find which is more realistic. 

To convert the dust optical depth gi ven by eq. 01 into an 
observable flux density, we follow lWvatti 12005) and assume the 
emission is from grains that behave as black bodies. Though 
the dust can in some cases be small enough to be hotter than 
a black body at a given distance from the star, tests using real¬ 
istic grain properties find that this approximation is reasonable 
(and if anything an underestimate of the flux) when calculating 
mid-IR emission; for A-type stars the smallest grains are large 
enough to behave approximately as blackbodies, and for later 
type stars the inefficient emission at wavelengths longer than 
the grain size is offset by their increased temperatures and con¬ 
sequently greater fluxes. 


2.2 Observed Kuiper belt analogues 

To connect observations of Kuiper belt analogues to the par¬ 
ent belt properties in the model, we must convert them to our 
model parameters ro and tq. The radius can be easily derived 
from the dust temperature of an observed disk assuming black- 
body properties, with the caveat that these estimates are typi¬ 
cally a factor of a fe w too low, depending on the spectral type of 
the star (e.g.lRodriguez & Zuckermanll20121 ; iBooth et aJ]|20 1 3l : 

IPawellek et al.1120141) . In any case, as we show below, an un¬ 

certain radius has a relatively small impact on the model pre¬ 
dictions. The optical depth can be estimated from the fractional 
luminosity f = Ldisk/T* = cr/(47rro), where a is the surface 
area of emitting dust, with the assumption of some fractional 
disk width via a = 27rroArro. Again, as we show below the 
habitable zone dust levels are relatively insensitive to the optical 
depth of the parent belt, so here we assume Ar = 0.5ro, and 
thus ro = 4/. 


For a sample of known Kuiper belt analogues we use stars 
from the Unbiased Nearby Stars (UNS) sample JPhi 11 ins et alj 
l2010l) . Most of these have been observed with Spitzer and/or 


Herschel (e.g. Rieke et al. 

2005J;|Trilling et al. 200$; Eiroa et al. 

20 ll; Thureau et al. 2014 

), and here we use disk properties de- 


rived from this far-IR photometry to show a representative disk 
population. The stellar photosphe res have been modelled with 
PHOENIX AMES-Cond models dBrott & Hauschildj2005i) and 
the excess fluxes with blackbodies to derive the dust tempera¬ 
tures and luminosities, whi ch are the quantities of interest here 
(e.g. iKennedv et al]l2012bH 3; iThureau et alJl2014h . These belts 
have radii of a few au to a few hundred au. and fractional lumi¬ 
nosities between 10 -6 and 10 -3 . 


3 HABITABLE ZONE DUST LEVELS 
3.1 An example 

Fig. U| shows radial optical depth profiles for a Sun-like star for 
both collision rates and several different parent belt parameters. 
The outer belt parameters are chosen to represent a relatively 
bright Kuiper belt analogue, and one near the current detection 



Radius (au) Radius (au) 


Figure 1. Examples of radial optical depth profiles interior to parent 
belts around a Sun-like star. The left-hand panel shows a linear radial 
scale while the right-hand panel is logarithmic. Four different parent 
belts are shown as dots, with ro = 30 and 90 au, and to = 5 X 10“ 6 
and 5 X 10 -4 . The darker lines show the high collision rate (k = 1/7) 
and lighter lines show the low collision rate (k = 1). The Solar Sys¬ 
tem’s Zodiacal cloud level (T zodi’) is also shown from 0.5 to 2 au. 
The dust optical depth in the habitable zone is near-constant, and de¬ 
pends strongly on the collision rate and weakly on the optical depth in 
the parent belt. 


limit. The optical depth of the Solar Zodiacal cloud (i.e. 1 ‘zodi’ 
unit of dust) is shown from 0.5 to 2 au dKelsall et al.ll 19981) . 

Dust is created in the parent belt with optical depth to, and 
then moves towards the star due to P-R drag, being depleted by 
collisions on the way. The collision rate depends on the optical 
depth, so eventually P-R drag dominates. The P-R drag time is 
oc 1/r so the optical depth approaches a constant level with 
radius. Thus, the optical depth interior to very tenuous parent 
belts (where collisions are unimportant) is also constant. 

For denser parent belts the interior levels become insensi¬ 
tive to the parameters. More dust drifts inwards from denser 
belts, but this extra dust is also collisionally depleted more 
rapidly. For fixed stellar parameters and parent belt radius, the 
interior optical depth reaches a maximum optical depth near the 
star a s tq is increased (i.e. no 1 and r(r) is independent 
of To. iWvattll2005l) . This insensitivity to parent belt properties 
is an important aspect of the model that is not expected in other 
scenarios, where stochastic behaviour and/or strong dependence 
on system architecture are expected to lead t o a wid e vari- 
ety of habitable zone dust levels (e.g. [Kennedy & Wvatt]|20 1 3l : 
iBonsor & Wyatt 2012). However, the inner dust levels depend 
strongly on the collision rate assumed (i.e. k). For the exam¬ 
ples in Fig.Q]the high collision rate predicts dust levels that are 
approximately 10 times the Solar Zodiacal level, and the low 
collision rate levels about 50 times larger. In either case, in the 
absence of effects that remove the inward moving dust, the re¬ 
gions interior to Kuiper belt analogues seen around other stars 
(with / > 10 -6 ) will have dust levels that are significantly 
higher than the Solar Zodiacal level. 

In the outer Solar System the dust level interior to the 
Edgeworth-Kuiper belt, with r ~ 10“ 7 dVitense et al.ll2012l ). 
is thought to be constant in to about 10 au, where most of the 
particles are ejected by Saturn and Jupiter (estimates of the frac¬ 
tion reaching Earth’s vicinity depend on grain size, and vary 
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from 0-20%. lLiou et alJll996l ; lKuchner & Starkll20 1(j|). The od - 
tical depth of the Asteroid belt is similar l lDermott et alJl2003 1. 
but most of the dust that spirals in towards the Earth is not lost. 
Given the uncertainties in their inferred optical depths and mod¬ 
els, both have been considered a sou rce for our Zodiacal cloud 
dLiou et alJl996l ; lGrogan et al J200 ll) . In reality, both contribute 
at some level, and material from Jupiter-family comets serves as 
a third and equally important source (i.e. dynamical r ather then 
P-R drag delivery fr om the Kuiper belt, Wh ipple et al.l 1 19671 : 
INesvomv et alJl2010h . Thus, external observations of the Solar 
System’s dust complement (that might not detect the gap be¬ 
tween ~3-5 au), might reasonably conclude that the Zodiacal 
dust originates in the Kuiper belt and is delivered by P-R drag. 
However, our knowledge of the actual architecture shows that 
planets can strongly influence the habitable zone dust level, both 
depleting it relative to simple expectations by removing dust as 
it spirals in, and enhancing it by allowing other delivery mech¬ 
anisms. 


3.2 Model predictions 

Fig- El shows predicted dust levels in the habitable zone for an 
A-type and a Sun-like star. We use the concept of the “Earth- 
equivalent insolation distance” (EEIDj as a simple definition of 
the habitable zone location, which is simply e qual to ^/L-^/Lp, 
in units of au (i.e. 1 au for the Sun. lRoberge et al.l2012l) . In each 
panel contours show the dust level at the EEID, relative_to the 
S olar System optical depth at 1 au of 7.12 x 10“ 8 (iKelsall et al.l 
11998b , for a range of parent belt optical depths and radii. That is, 
at the EEID (but not necessarily elsewhere) the se contours are 
the same as the “zodi” units defined in lKennedv et al.1 120151 ) for 
the purposes of LBTI modelling. Sub-panels for both the high 
and low collision rate models are shown. The symbols mark 
known Kuiper bel t analogues seen around nearby stars from the 
sample of dPhillips et alj|201(]t ). including a ‘tail’ that reflects 
the approximate systematic uncertainty in the dust radial loca¬ 
tion. Vertical uncertainties are relatively small (roughly a factor 
±2) as any model that fits the disk photometry will have a simi¬ 
lar integrated luminosity. A representative detection limit for the 
parent belts around nearby stars with mid and far-IR photometry 
(24 and 100 )im) is shown (computed as in lWvattl2008l) . where 
detections are possible above the dashed line, and differs for A 
and G-type stars because their luminosities are different. At low 
parent belt densities the dust level contours are roughly constant 
with radius because P-R drag dominates. At high densities col¬ 
lisions dominate and the habitable zone dust level depends only 
on the radial proximity of the parent belt to the habitable zone. 

Thus, P-R drag of dust from known Kuiper belt analogues 
can result in significant HZ dust levels relative to the Solar Sys¬ 
tem level, regardless of the assumed collision rate. The results 
clearly depend on the adopted collision rate; the low collision 
rate implied by the KIN results suggests that the HZ dust level 
is typically 50-500 times the Solar Zodiacal level for detectable 
parent belts around Sun-like and A-type stars. The higher colli¬ 
sion rate results in typical levels closer to 10-100 times the Solar 
Zodiacal level. 

The nature of the photometric detection limits means that 
there is a region of the parameter space in all cases where HZ 
dust levels in excess of a few tens of times the Solar System 
level can arise from belts that cannot be detected by photom¬ 
etry. These parent belts tend to be close to the habitable zone 


however, and therefore decay by collisions faster than belts that 
are farther out, mostly due to higher relative velocities. The ap¬ 
proximate maximum optical depth of a parent belt in collisional 
equilibrium is shown in each panel, for average stellar ages of 
500 Myr and 5 Gyr, using plan etesimal properties that fit their 
obser ved brightness evolution jWvatt et al.ll2007bl ; iKains et akl 
l201lh . The vertical position of these lines varies as l/f ag e, and 
for a given age few belts are seen or expected above such lines 
(e.g. I Wyatt et al.112007ah . Thus, the space where belts are al¬ 
lowed by collisional evolution but are not detectable and result 
in significant HZ dust levels is relatively small, and decreases in 
size as the age of the star increases. 


4 DISCUSSION 

Poynting-Robertson drag is one of many ways to produce habit¬ 
able zone dust, the main alternatives being material produced lo¬ 
cally in collisions or dynamical comet delivery by planet scatter¬ 
ing. Thus, where a Kuiper belt analogue is known the exo-zodi 
levels discussed above are maximum P-R drag induced dust lev¬ 
els that can be compared with observations to draw conclusions 
about the dust origin. That is, dust levels significantly greater 
than predicted require a different origin. 

For example, the most extreme system known to harbour 
warm and cool dust populations is the F2V star p Cor vi (e.g. 
IWvatt et al J 20051 : 1 Smith et all2009l : lDefrere et all2015h . where 
the location and optical depth of the warm belt is roughly 1 au 
and 10~ 3 , and 150 au and 10” 4 for the cool belt. Fig. |T] shows 
that the HZ optical depth expected from any parent belt from P- 
R drag is lower than 10 -5 , and therefore that the warm dust is 
not delivered from the outer belt by P-R drag. Indeed, the only 
systems where such bright warm dust may originate in cool belts 
are those with low stellar luminosities, where collisions during 
inspiral do not result in removal, and thus the interior optical 
depth re mains approximately the same as in the parent belt (e.g. 
e Eri. lReidemeister et alj20lH ). In general therefore, detections 
of warm dust well above the predictions of the P-R drag models 
around Solar and earlier-type stars require a different origin. In 
the case of t) Crv the favoured scenario is that comets are be- 
ing scattered inward from the cool b elt by unseen planets (e.g. 
IWvatt et all2007al : lLisse et all2012h . 


4.1 Habitable zone dust levels 


It seems probable that all stars host cool outer debris disks, with 
only the bright est 20-25% being detectable with current meth¬ 
ods (e.g. Trilling_etalJ 20081; Su_gtaL 20061 ; lEiroa et alj|2013l ; 
ISierchio et al.ll2014l ; lThureau et alj|20l4l) . Indeed, this was the 
assumption made in various population models th at rea son ably 
reproduce the statistics of such disk s (e.g. IWvatt et al.ll20 07b; 
Kains e t al.l201 HlGaspar et alfeoHii) . Given the inevitability of 
P-R drag, this population of cool outer belts implies the exis¬ 
tence of a corresponding population of warm dust around all 
stars. Combined with the results in Fig. [2] the 20-25% of stars 
with detected cool belts can therefore be used to make some 
useful statements about HZ dust levels. 

For A-type stars. Fig. [2] shows that detectable cool belts 
result in HZ dust levels at least 20-100 times greater than in 
the Solar System (i.e. most known outer belts are above these 
contours). Thus, among the remaining 75-80% of stars the HZ 
dust levels due to P-R drag are below these levels. For G-type 
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Figure 2. Dust levels at the EEID for an A-type (3.7 au) and a Sun-like (1 au) star for a range of ro and to. Contours show the optical depth relative 
to the Solar System level at 1 au of 7.12 X 10 —8 . Dots mark known Kuiper belt analogues around nearby stars and connected lines show the probable 
increase in size due to non-blackbody emission from small dust. An approximate photometric detection limit is shown by the dashed line. Dot dashed 
lines show the maximum fractional luminosity as a function of radius for stellar ages of 500 Myr and 5 Gyr (as labelled). For each panel, the left 
sub-panel shows results for the low collision rate (k = 1), and the right sub-panel shows results for the high collision rate (k = 1/7). The left end of 
the radial scale stops at the EEID in each case. 


stars the predicted HZ dust levels are lower, and the 75-80% 
or so of Sun-like stars with no known outer belts have HZ dust 
below 10-50 times the Solar System level. 

Various studies have estimated the impact of HZ dust on 
future missions that will attempt to detect and characterise 
Earth-like planets around other stars. The two most recent stud¬ 
ies seem to reach a consensus that, with the appropriate survey 
strategy and under the assumption of smooth exo-zodi, this im¬ 
pact is a relatively weak function of dust brightness, and that for 
levels below a few tens of times the Solar System level a mis¬ 
sion will not be seriously hindered dStark et al.l 120141 : iBrownl 
120151) . The possible impact of non-axisymmetric structure is 
less certain however, since inhomogeneities in exo-Zodi surface 
brigh t ness may themselv es be mistaken for planets (e.g. lLavl 
120041 : IDefrere et alj|2010f ). Such structures may of course also 
be sign po sts of th e dynamical influence of unseen planets (e.g. 
IStark & Kuchneiibopsl ). Whether such structure presents a ma¬ 
jor barrier depends both on the typical level of non-axisymmetry 
and the specifics of the Earth-imaging mission, but could result 
in a tolerable exo-Zodi level closer to ten times the Solar System 
level. 

Based on our HZ dust predictions, Earth-imaging efforts 
for Sun-like stars are unlikely to be seriously affected by dust 
coming in from undetected outer belts, but most of the 20-25% 
of systems with detected outer belts would be considered unsuit¬ 
able targets for this reason. For A-type stars such a statement 
is less certain, as it is possible that 10-30 au belts just below 
the level of detectability (whose frequency is unknown), will 
contribute to HZ dust levels. Collisional depletion of the parent 
belts over time means that this issue can be mitigated to some 
degree by avoiding young stars. It is unlikely that such belts 
make up a large fraction of the population however, so the frac¬ 
tion of unsuitable systems is probably not significantly greater 
than 25%. 

Thus, while a non-negligible fraction of stars are expected 
to be poorly suited to Earth-imaging/characterisation due to P-R 
drag from outer belts, the detection limits of extant far-IR sur¬ 


veys means that among stars within a few tens of parsec most 
of these systems are already known. As noted above however, 
P-R drag is by no means the only origin of habitable zone dust. 
Thus, in the absence of other effects (see next section) the lev¬ 
els predicted by detections of cool outer belts represent a lower 
limit on dust in a given system. 


4.2 Inference of unseen planets 

The inward transport of dust from a known parent belt by P- 
R drag is inevitable, so additional processes must also be in¬ 
voked if HZ dust levels could be shown to be significantly 
lower than expected. That is, something may be needed to re¬ 
move the dust somewhere between the parent belt and where 
detection was attempted. By analogy with the Solar System, 
the most obvious reason is of course accretion or ejection of 
that dust by intervening planets i Liou et all 1 9961 ; iLiou & Zookl 
1 1 999l : [Wyatt et al.ll 1 999l : fvitense et aTll20 1 2t) . Whether the dust 
is accreted or ejected will to first order depend on the ratio of the 
escape velocity from the planet to the local escape velocity from 
the star. Thus, a further implication is that close-in planets (e.g. 
those that transit their stars) are also those most likely to accrete 
dust, rather than eject it, perhaps with observable effects. The 
limits set by the KIN were unfortunately not stringent enough 
to require the explanation of dust removal by other processes 
for any systems, their Fig. 10 shows that the non-detections are 
close to or above the levels predicted by the low collision rate 
model. 

Perhaps the most promising system where such a test could 
be made is around the “retired” A-type star k CrB. Though the 
spectral type is K0, the star is a sub-giant and for our purposes 
can be treated as an A-type star because stellar luminosity and 
mass are the important model parameters here. This system has 
a bright outer dust belt that extends from 20-40 au out to 165- 
220 au depending on the disk model dBonsor et al .. 2 0131) . and at 
least two companions. One at 2.7 au has been well characterised 
by radial velocity Ijohnson et al.l2008l) . and another is inferred 
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to exist at greater distance due to the existence of a residual ac¬ 
celeration in the radial velocity data dBonsor et al.ll2013l) . With 
only a constant acceleration seen over 8 years, the outer com¬ 
panion lies beyond about 7 au, and therefore beyond the EEID 
of 3.5 au. This companion is inferred to be more massive than 
Jupiter, and resides at a radial distance where the escape veloc¬ 
ity from the star is lower, and can therefore easily eject particles 
that interact with it. Thus, if this companion lies interior to the 
outer dust belt it will eject much of the dust that would other¬ 
wise reach the habitable zone. Given confidence in the P-R drag 
model, a sufficiently constraining non-detection of warm dust 
around k CrB could therefore be used as further evidence of the 
second companion’s existence. 

In the Solar System, the Asteroid belt and Jupiter family 
comets provide alternative habitable zone dust sources, and it 
seems largely a coincidence that the Zodiacal cloud has ap¬ 
proximately the level expected for P-R drag from Kuiper belt 
dust. As noted in section [TT] an external observer with similar 
but moderately more sensitive instruments than ours might con¬ 
clude that the dust in the inner Solar System is consistent with 
that expected from P-R drag from the Kuiper belt, and therefore 
that invoking unseen planets is not necessary. Thus, while ad¬ 
ditional processes and dust source regions may cause existing 
planets to be missed, they would not cause non-existent planets 
to be inferred. 


4.3 Detectability with LBTI 


An uncertainty in the above predictions for HZ dust levels is 
which of our two models is more correct. That is, a significant 
step towards improving the model would be an empirical cali¬ 
bration of the parameter k across a range of stellar spectral types 
and parent belt properties. We now show that these models will 
be strongly tested by the LBTI, and hence aid further model de¬ 
velopment. The LBTI will survey a sample of ~50 nearby stars, 
with a wide vari et y of spe ctral types, to look for habitable zone 
dust ^Weinberger et al.ll2015h . The sample includes stars known 
to host cool outer dust belts; a few of these are Sun-like stars 
but most are earlier A-types. 

Briefly, the LBTI is a nulling interferometer that operates 
in the mid-IR (11 /im), so interferes the light collected by each 
of the two LBT mirrors to attenuate the starlight. Off-axis mid- 
IR emission at an angular scale of greater than ~40 mas, for 
example from warm dust, is transmitted. The attenuation of the 
starlight avoids the problem of disentangling the dust emission 
from the stellar emission, which is the limiting factor for dis¬ 
covery of habitable zone dust using purely photometric meth¬ 
ods (i.e. sets the dashed lines in Fig. [2). In practice, LBTI ob¬ 
servations involve various calibration steps, and for a descrip¬ 
tion of these we refer the reader to iDefrere et ali J2015I) . The 
key point is that an LBTI measurement, known as the “cal¬ 
ibrated null depth” or “null excess”, is the ratio of the disk 
flux that is transmitted though the LBTI fringe pattern to the 
star+disk flux. Thus, because the disk is almost always much 
fainter than the star, the null excess can be thought of as anal¬ 
ogous to the disk to star flux ratio, but with the disk flux at- 
tenu ated by the sky projection of the LBTI trans mission pat¬ 
tern (i Mi 11 an-Gabet ~et al fed 1 ll; Defrere et al.ll2015l) . We use the 
model outlined by iKennedv et alj G015l l. and assume face-on 
disks to convert our model dust levels to null excesses. 

Null excesses are shown as contours in Fig. [3] again for 


a range of parent belt properties, for both Sun-like and A-type 
stars and the low and high collision rates. Again, known outer 
belts around nearby stars and photometric sensitivity limits are 
also included. The projected LBTI 3 a null excess sensitivity is 
3 x 10~ 4 , so these contours, and a higher level of 10~ 3 , are 
shown. Because the null excess depends on the angular scale 
of the disk, the contours show pairs at two different distances 
(TO and 30 pc). The null excesses always increase with opti¬ 
cal depth, so for each distance the contours show the minimum 
parent belt optical depth for warm dust to be detectable at that 
level. 

Looking first at the low collision rate model (left panels), 
all of the known parent belts lie in regions where the warm dust 
is detectable by the LBTI. In addition, for A-stars there is a large 
swath of parameter space for LBTI detection of warm dust from 
parent belts that are otherwise undetectable by photometry. A 
similar region exists for Sun-like stars, but covers less parent 
belt parameter space. For this model warm dust is therefore pre¬ 
dicted to be detectable with LBTI, easily in most cases, for all 
systems with known Kuiper belt analogues. Thus, if the KIN 
results of a direct connection between warm and cool dust is 
the result of P-R drag and the low collision rate model applies 
across AFG spectral types, the LBTI will strongly confirm this 
result. 

In cases where no parent belt has been seen, warm dust 
originating in such a belt could still be detected. The origin 
of this low level HZ dust will be ambiguous; mere detection 
will be sufficiently difficult that further constraints on the ra- 
dial extent, which could provide clues to the origin, will be poor 
dKennedv et al.ll20l5l ; lDefrere et ro 

For the high collision rate the results are very similar for 
A-type stars; nearly all systems with known parent belts should 
result in LBTI detections. Detections of dust from undetected 
parent belts also remains possible. For Sun-like stars however, 
only systems with parent belts with radii <15 au should result 
in LBTI detections. 


5 SUMMARY AND CONCLUSIONS 

The development of new instruments means that P-R drag can 
no longer be considered an insignificant effect for transporting 
dust to regions interior to the Kuiper belt analogues seen around 
other stars. While the point that the contrast between a parent 
belt and the inner regions is high remains, the absolute dust 
levels interior to known parent belts can be hundreds of times 
higher than seen in the Solar Zodiacal cloud. We have explored 
these levels using two flavours of an analytic P-R drag model, 
connecting the properties of known Kuiper belt analogues to ex¬ 
pected habitable zone dust levels. 

Habitable zone dust levels greater than a few tens of times 
the Solar System level will be detrimental to future missions to 
discover and characterise Earth-like planets around other stars. 
Currently detectable Kuiper belt analogues are predicted to re¬ 
sult in HZ dust levels similar to this limit. Thus, mid and far-IR 
debris disk surveys have already identified most of the 20-25% 
of stars where P-R drag from outer belts could seriously im¬ 
pact Earth-imaging. The existence of other processes that gen¬ 
erate/deliver HZ dust of course means that the 20-25% is a lower 
limit. 

A caveat to such predictions is that intervening planets may 
prevent this dust from reaching the habitable zone, instead being 
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Figure 3. LBTI null excesses for a range of parent belt parameters, for two different sensitivities (3 X 10 —4 and 10“ 3 ), and two stellar distances (10 
and 30pc, contours as labelled). The top row shows models for a Sun-like star, and the bottom row shows models for an A-type star. Left panels show 
models for the low collision rate, and right panels show models for the high collision rate. Dots mark known Kuiper belt analogues around nearby 
stars and connected lines show the probable increase in size due to non-blackbody emission from small dust. Dashed lines show approximate detection 
limits for parent belts around nearby stars. The null excesses always increase with parent belt optical depth, so for each distance the contours show the 
minimum optical depth for warm dust to be detectable at that level. 


ejected or accreted by intervening planets. Thus, in systems with 
known outer belts the non-detection of such warm dust could, 
given sufficient confidence in the models, be used to infer the 
existence of unseen planets. Currently the models are uncertain 
however, so the goal of future work should therefore be first 
to empirically calibrate the P-R drag models (e.g. with LBTI 
observations) to aid model development, and then to interpret 
the observations within the framework of such a model. 
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